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An improved formula is derived for accurately computing the near-zone magnetic field of a small 
circular transmitting loop antenna. Such a field can serve as a reference standard for calibrating 
field-strength meters employing small receiving loop antennas in the frequency range 30 Hz to 30 MHz, 

This formula includes correction terms for frequency (due to the finite time of propagation), as 
well as corrections for the finite radii of both the transmitting and receiving loops. Other formulas 
appearing in the literature often fail to include such corrections which can result in errors of up to 
20 percent and more in computing standard-field values. 

The NBS formula is derived by expanding the integrand of the retarded vector potential into an 
infinite series of splierical Hankel runctions of increasing order. Th( resuking series expression is in 
error by less than 0.2 percent, is rapidly converging and simple to use without recourse to a table of 
functions or a computer. 
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1. Introduction 

The purpose of this paper is to present a derivation of 
an improved high-frequency formula intended for use 
in accurately computing the near-zone magnetic com- 
ponent of the electromagnetic field of a small circular 
transmitting loop. Such a field can serve as a reference 
standard for calibrating field-strength meters employ- 
ing small receiving loop antennas in the frequency 
range 30 Hz to 30 MHz. The new formula includes 
correction terms for frequency (due to the finite time 
of propagation), as well as corrections for the finite 
radii of both the transmitting and receiving loops. 
Formulas used by other workers in the field often do 
not include such corrections. This new formula is a 
rapidly converging infinite series. It is simple in form 
and easy to use by workers in the field without recourse 
to a table of functions, or to a computer. 

An earlier formula, which prompted this derivation, 
has been used during the past 15 or 20 years at the 
National Bureau of Standards [1].^ It has served as the 
basis for the NBS calibration service for commercial 
and military CW field-strength meters over the above 
frequency range. This formula was a single-term ex- 
pression and was first obtained from a detailed numeri- 
cal analysis of the existing low-frequency formulas, of 
which there were many in the literature [2]. Most of 
these were in the form of an infinite series. Many were 
slowly converging and required the use of at least two 
or three terms of the series to yield the accuracy 
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achieved in the single term of the earlier NBS formula, 
over the usual range of the parameters involved. 

From the form and symmetry of this single-term 
expression, it became apparent that it might be the 
leading term of a rapidly converging infinite-series 
expression for the magnetic field of the loop. Efforts 
were then directed toward finding an analytical solu- 
tion in series form having the earlier formula as its 
leading term. A number of methods were tried without 
success. The one that finally achieved the desired re- 
sult is based on expanding the integrand of the re- 
tarded-vector-potential into an infinite series of spheri- 
cal Hankel functions of increasing order. While this 
approach must have occurred to other workers in the 
past, the author has been unable to find any record of a 
high-frequency derivation in this particular series 
form in the literature. 

All time-varying quantities used in this paper are 
expressed in terms of their rms values. Rationalized 
mks units are used throughout. 

2. Calibration and Use of a Small Loop 

The physical arrangement used at NBS for cali- 
brating small receiving loop antennas is shown in 
figure 1. The standard magnetic field is produced by 
a small, single-turn, circular transmitting loop of 
radius, r\ meters, in which a current, / amperes, is 
flowing. The small receiving loop being calibrated 
has a radius, r-y meters, and is positioned coaxially 
with respect to the transmitting loop. The two loops 
are spaced an axial distance, d meters, apart during 
the caHbration. It is assumed that the distance of 
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Figure 1. Coaxial arrangement of transmitting loop antenna 
{radius Y\ meters), and receiving loop antenna {radius Xi meters), 
spacing d meters, in the NBS calibration setup. 

An RP" current, / amperes, flows in the transniittinji loop. 



separation, d, is small compared to the wavelength, 
X, and that near-field conditions prevail. 

The response or indication of the field-strength 
meter is directly proportional to the voltage induced 
in the receiving loop. The induced voltage is, in turn, 
proportional to the time-rate of change of the total 
flux cutting the loop. This is given by Faraday's law 
of electromagnetic induction, written here for the 
sinusoidal time-varying case (with the time factor 
e^^^ understood). 



E ' ds2 = - 



-J^ 



// 



B • dS-> 



(1) 



The voltage induced in the loop is given by the left- 
hand side of eq (1), which is the line integral of the 
tangential component of the incident electric field, E^ 
around the periphery, so, of the receiving loop; B 
is the magnetic-flux density at any point on the plane 
surface, ^2, bounded by the circular receiving loop; 
CO = 277 f, and j= V^. 

When making a field-strength measurement, one 
does not know the spatial distribution of B over the 
area, ^2, and it cannot be resolved by the loop. The 
response can therefore only be expressed in terms of 
the average field over the area of the loop, and eq (1) 
becomes 



V=-jcoB, 



S-z, 



(2) 



where V is the jriduced voltage and B^^ is the normal 
component of B averaged over the area, ^2, of the 
receiving loop. This assumes that the relative orienta- 
tion of the receiving loop with respect to the electric- 
field vector, E, is such that the electric-dipole response 
is negligible [3]. 

3. Retarded Vector Potential of a Circular 
Transmitting Loop 

The basic far-field properties of a small transmitting 
loop antenna are usually derived for conditions under 
which second-order corrections to the field due to the 
finite radii, r\ and ro, of the transmitting and receiving 
loops, can be neglected. In nearfield use, however, 
Ti and r-z are usually not negligible compared to the 



distance of separation, d, between the loops. There- 
fore, corresponding corrections must be provided in 
the field magnitudes involved by including both of 
these radii in the derivation. 

The expression for the axial magnetic component 
of the near field of the electrically-small (277/ 1 < k) 
transmitting loop antenna will now be derived. It will 
be assumed for the present, that the loop current is 
constant in amplitude and phase around the loop, and 
that the radius of the wire with which the loop is wound 
is neghgible compared to the loop radius, n, itself. 
The effect of a non-uniform current distribution on the 
accuracy of the resulting field formulas at frequencies 
above 5 or 10 MHz will be discussed later in section 5. 
The effect of the finite time of propagation will be 
included in the derivation. _ 

The magnetic-flux density, B, at any point, P, on 
the surface, ^2, of the receiving loop, is given by the 
curl of the retarded vector potential, ^, at the same 
point [4], i.e.. 



B = VxA. 



(3) 



The geometry involved is shown in figure 2. The aver- 
age value of the normal component of 5 will be given by 



B. 



where [4] 



-ill 



VxA ■ dSi, 



(4) 



R ^^" 


(5) 


pace {fi = 4^7T ' 10"' 


, henry/ 



meter) 

/= transmitting-loop current, rms amperes 
/3 = wavelength constant (P = 2ttI\) 
X = free-space wavelength, meters 
/= frequency, hertz 



It can easily be shown that the vector potential will 
have only a component, A^, which is azimuthal at the 
point, P. Let P be chosen arbitrarily to lie in the xz 
plane where (^ = 0. Equidistant line elements of length 
ds\ at +(/) and —cj) are paired. These elements can 
each be resolved into a component, ds\ sin (/>, parallel 
to the xz plane and a component, dsi cos (/), normal to 
the xz plane. In taking the resultant, the former com- 
ponents will cancel, while the latter (azimuthal with 
respect to the point P) will add. Equation (5) can 
therefore be rewritten [5], letting dsi = r\dcf). 



A, 



jUr, 

' 277 



I R 



cos (/) d4>. 



(6) 



320 



By Stokes's Theorem f6|, the integral of the normal 
component of the curl of T over the surface ^2 bounded 
by the receiving loop can be transformed mto the line 
integral of the tangential (component of .4 around its 
periphery, i.e.. 



// 



VxA ' (IS>= (p A4F)ds2, 



(7) 



where A^t>{P) is the value of the vector potential at the 
point, P, now assumed to He on the periphery of the 
receiving loop of radius, p = r2. Substituting eq (7) 
in eq (4) gives 



B. 



1 



A^{P)ds2. 



(8) 



Since the receiving and transmitting loops are coaxial, 
it can be seen from the circular symmetry involved, 
that Afi,(P) will be constant around the periphery of 
the receiving loop, so that eq (8) becomes simply 



B,, = ^A4P). 



(9) 



Substituting eq (6) in eq (9), letting B ,,,, = flH ,,^. and 
^2 = TTrs, gives 



//av = 



In 

77T: 



1 1 ~ii~ 



cos (f) dcj). 



where now 



R=[d'-^n-^4- 2r, r. cos 1 '/^ 



(10) 



(11) 



as can be determined from figure 2. H.^^ is the average 
normal component of magnetic field strength over the 
area of the circular receiving loop. 

4. Evaluation of the Complex Integral 

The problem, now, is to evaluate the integral in 
eq (10). Maxwell [7] was the first to accomplish this for 
the static case (^==0), using elHptic integrals which 
can be expressed in several series forms [8]. Numerous 
other workers have since evaluated this integral in 
various ways for the static case. Literally, over one- 
hundred different series expressions appear in the 
literature [9] for the mutual inductance between 
coaxial circular filaments, to which the present prob- 
lem is directly related. 

To the above will now be added still another series 
solution, this time also including the effect of retarded 
time. Equation (10) will be evaluated using the rela- 
tionship [10] 



jfi/i 



-m 



h'vm\ 



(12) 



where h^l\f3R) is the zero-order spherical Hankel 
function of the second kind. Such a function is often 
used to represent an outward-traveling ( + ) wave. 
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Fl(;iiRE 2. Geometry used for determining the vector potentidl ami 
near-zone magnetic field of a small circular transmitting loop. 

Tlic l(>(i|) lies in a plane normal to the xz plane with its center at the orijiin of coordinates. 



whereas a similar function of the first kind, /?\V(/3/?), 
would have been used to represent an inward-travel- 
ing (— ) wave. The function, h^'^i^R), will now be 
expanded into an infinite series of higher order spheri- 
cal Hankel functions. These will be substituted back 
into eq (10) and integrated term by term. 

The expansion can be accomplished using the fol- 
lowing relationship based on Erdelyi |11|, 



/i',V(y3/?) = A<?,'(V;7+^) = 2 



I 



A',t/(V«). 



(13) 



where m is an integer. This series converges when 
|u|>|i^|. Let 





Vu + 1; = y3 (/?g - 2r, r2 cos </) ) "2 , 


(14) 


where 


Rl=cl} + fi + ri. 


(15) 


Let 


u = ^mi. 


(16) 


and 


v = — ft^{2r,r2 cos </)). 


(17) 



Substituting eq (16) and eq (17) in eq (13) gives 
1 



h^iKm=^ 






/il2)(^/?o) cos- 0. (18) 



Substituting eq (18) back into eq (10) gives 



(19) 
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The infinite series in eq (19) will converge [11] when 
/?§ > 2rir2 if (/) = or if = 77. Furthermore it can be 
shown [12] using, for example, the Weierstrass M 
Test, that eq (19) will converge uniformly in the inter- 
val ^</)^7r. Equation (19) can then be integrated 
term by term, yielding 



//. 









Jo 



5'"+i cl)d(t), 
(20) 



The integral in eq (20) can be evaluated from the follow- 
ing relations [13]: 



jo 

/' 

jo 



cos^^ (f)d(j) = 0, 



cos^^ (fxlcf) - 



1 • 3 • 5 



2 •4-6 



-TT, 



(if p is 


odd) 




(21) 


(if J) is 


even) 




(22) 



When these results are substituted into eq (20), with an 
appropriate change in the index of summation, eq (20) 
becomes 



//av = 






1 



1 • 3 • 5 



'^ (2m+l)!2-4-6 



2m + 1 



2m + 2 



Ro 



^,.At^^o). 



(23) 



The spherical Hankel functions in eq (23) can now be 
changed into their exponential polynomial forms as 
given by eq (All) to eq (A16) in the Appendix. When 
this is done, eq (23) can be written 



"-=^M' +-''^«»] +?(1)' \l+jl3R.-^^p-^Ri-jl-^l5mS 



+ 



315/rir2 
64 [rI 



1 ^mo - 13 fsm,^ -J If fS^S + ^ f^^R^ + J ^ f^'RS 



+ 



-Jfif^o 



amperes/meter. 



(24) 



Under the conditions that jS/?o ^ 1.0, and rir-zIRl 
^ 1/16, terms in fBRo higher than the first power can 
be ignored without introducing an error in eq (24) of 
more than 0.2 percent. Therefore, if the term (l+//S/?()) 
is factored out of eq (24), the magnitude of eq (24) can 
be written as 



\H.J 



' 2ttR§ 



l + '-^M 



+ 



315 
64 






d + fi-'Ro') 



(25) 



This is the desired expression for the normal com- 
ponent of the magnetic field averaged over the area, 
S2, of the receiving loop. 

The effect of the neglected terms in l3Ro on the 
accuracy of eq (25) decreases rapidly as l3Ro is pro- 
gressively decreased below 1.0, other factors remaining 
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the same. At the lower frequencies where correc- 
tions due to fSRo can be neglected entirely {^Rq < 1), 
the mathematical accuracy of eq (25) is limited only 
by the number of terms used in the infinite series. 
For the conditions previously jiivcn, the first cor- 
rection term in the infinite series of eq (25) will con- 
tribute less than 1 percent, so that for this case eq (25) 
can be further simplified to 

|/yavl=2^^^ ( 1 +f^~R;iy'\ umperes/meter. (26) 

It should be noted that eq (26) is essentially identical 
to the expression for the axial component of the 
ma^inetic field of a classical, infinitesimal magnetic 
dipole. The difference is that the distance factor, /?(), 
in eq (26) includes correction terms for the finite radii, 
Ti and r-z, of the transmitting and receiving loops, 
respectively, as given by eq (15). 

As previously stated, the response of a field-str(Mi<i;tli 
meter ein|)l()ying a small-loop antenna is directly 
pro|)orti()nal to the avtMage normal component of the 
magnetic fi(dd strength, H.^y, incident on the loop. 
However, ior many years it has been customary to 
express the calibration oi such an instrument in terms 
of the ('(juivalent electric field strength, A\ that would 
be associattul with H for the case of a uniform plane 
wave (wher<^ E=\2i)7TH). \\us relationshij), wIhmi 
substituted in eq (26) with Si = 7rr^, gives 



l^avl 



6()7r/f/ 



(i+iS^/^,r 



i«/^ 



volts/meter. (27) 



Needless to say, this maimer of expressing the re- 
sponse will be valid o/?/v when the field-strength meter 
is used to measure a uniform plane wave. The meas- 
urement will be of questionable value for (lU cases in 
which the local |A///| vdUo departs from the free-space 
value [14]. 

Standard-fi(^ld lornuilas appearing in the literature 
[15] to [17] often fail to correct for the finite radii of 
the transmitting or receiving loop antennas, or for 
frequency. It is of interest to evaluate the errors 
resulting from the failure to apply these corrections, 
all of which are included in the formulas derived in 
this paper. 

Failure to correct for the radius, r-z, of the receiving 
loop will result in an error of from 1.0 percent (for 
d=l.25 m, r2 = 0.1 m), to 20 percent (when d=0.5 m, 
r2 = 0.2 m). Failure to simultaneously correct for the 
radius, n, of the transmitting looj) will approximately 
double these errors. 

The value of field strength given by eq (25) to 
eq (27) is essentially independent of frequency at 
the lower frequencies (where jBRo < 1). Failure to 
apply the frequency correction indicated, at the higher 
frequencies, will introduce an error of, e.g., 1.0 per- 
cent at 5 MHz, which increases to approximately 27 
percent at 30 MHz (when d=l.25m, ri = r2 = 0.1m). 
If the transmitting loop circumference exceeds 
ap|)roximatcly \/16, tlie effect (da nommifoiin current 




Figure 3. Geometry used to analyze the current distribution on a 
circular transmitting loop antenna. 

The loop is treated as a circular balanced transmission line, fed at points 1, 2. and short- 
circuited at the receiving end (point 3). Distance along the loop periphery, / = 6i^^ meters, 
is measured from point 3. 



distribution will also have to be taken into account, 
as discussed in the next section. 

5. Loop-Antenna Current Distribution 

In the preceding derivation, the loop current distri- 
bution was assumed, for simplicity, to be constant in 
amplitude and phase around the loop. Actually, a 
standing wave of current exists on the loop, and 
the assumption of a constant current will not be valid 
unless the loop is small electrically (i.e., unless 
277/-, <^ X). 

The complete solution to this problem for an 
arbitrary current distribution [18J is quite complex. 
However, the effect of a nonuniform current distribu- 
tion at frequencies from roughly 5 to 30 MHz can be 
determined from the following approximate analysis. 
The single-turn loop can be considered as a circular, 
balanced transmission line fed at points 1, 2, as 
indicated in figure 3, and short-circuited at its receiving 
end (point 3). 

If the radius of the wire with which the loop is wound 
is small compared to the loop radius, n, it can be 
assumed that the current distribution around the loop 
is, to a first approximation, of the hyperbolic-cosine 
form [19], 



/ = /() cosh yl 



(28) 



where 



/() = reference current at point 3 (fig. 3) 
y = propagation constant = a -\-JI3 
a = attenuation constant, nepers/meter 
/3 = wavelength constant = 27r/X, radians per meter 
X= free-space wavelength, meters 
/= peripheral distance, in meters, measured from 
point 3. 
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The loop will be considered sufficiently loss free 
that a can be assumed negligible compared to ^. 
This is a reasonable assumption for the type of loop 
used, so that eq (28) can be written 



I = h cosh {j^l) = /() cos 131. 



(29) 



As can be seen from eq (29) the current will have a 
constant phase around the loop for this case. 

For the size loop of interest here, the magnitude of 
the magnetic field strength given by eq (24) to eq (26) 
will be, to a first approximation, proportional to the 
average value of the current around the transmitting 
loop. The average current is given by 



-' av 






131(11 = 1, 



sin j8/i 



(30) 



where /i = half the loop circumference, meters and 
I = 6?\ meters from fi<>;ure 3. 



Table 1. Percent difference between lyv and lo 
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Figure 4. Relative current distribution around a 
transmitting loop antenna having a radius ri = 0. 
quency (=30 MHz. 

The current distribution, based on an analysis by T. T. Wu [20J. i 
assumed cosine distribution determined from e(i (29) in the text, where 
A wire radius a = 0.001 ni was used ui the Wu analysis. 
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I m, for a fre- 

s compared with an 
/ = Oil from fi<:ure 3. 



6. An Alternate Method for Evaluating the 
Integral 

The infinite series of eq (25) can be easily verified 
for tlie static case {jS = 0) using a much simpler method. 
This is done by rewriting eq (10) in the form 

//av = -^ r Fict>) cos (/)J(/), (31) 



Table 1, which is based on eq (30), shows the 
approximate percentage difference between I^y and 
/o for several sizes of low-loss loops {a < j8), with 
their circumferences expressed in fractions of a 
wavelength, lirrxlX. As can be seen, the correction 
due to the nonuniform current distribution will be 
less than 1.0 percent at the highest operating fre- 
quency, provided the loop circumference does not 
exceed \/16. This correction decreases rapidly 
with decreasing frequency, becoming less than 0.2 
percent at frequencies below 15 MHz for a loop 
having a radius ri = 0.1 m. 

These corrections, based on an assumed cosinus- 
oidal current distribution around the loop, are in quite 
close agreement with similar corrections based on a 
recent theoretical analysis by T. T. Wu [20]. The dif- 
ference in the correction values obtained by the two 
methods in less than 0.02 percent for the conditions 
involved here. It is interesting to note that for this case 
the loop current distribution obtained by Wu reduces 
essentially to the cosine form given by eq (29) above. 
This is shown in figure 4 for a loop having a radius 
ri = 0.1 m operating at a frequency /= 30 MHz. 



where 



FW 



2rir2 cos (/) 

^0 



(32) 



The binomial expansion of eq (32) is then substituted 
back into eq (31) and the result integrated term by 
term, yielding the infinite series of eq (25) directly, 
provided Kq > 2r\i^2 as before. The first few terms of 
the binomial expansion of eq (32) are as follows: 



F(<i>) 



l + l^l cos </) + - 



+- 



Rl 



5 /ri/V 



2\Rl 



COS- (/) 



2\Ri 



COS'^ (/)-h 



35 



rxi'i 
til 



COS^ (\> 






(33) 



Substituting eq (33) back into eq (31) and perform- 
ing the required integration, making use of eq (21) 
and eq (22), yields for the static case 
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//=, 



AS, 






(34) 



which is identical to the iiifinih- series in eq (25). 

7. Summary and Conclusions 

An improved hi<»;li-iV(Miueiiey (onnula has been de- 
rived for accurately eoniputin'i the near-zone magnetic 
field of a small (27rri < X), circular, transmitting 
loop antenna. Such a field is used as a reference 
standard at NBS for the calibration of field-strength 
meters [1], [14] employing small receiving loop an- 
tennas in the frequency range 30 Hz to 30 MHz. 

The new formula is in the form of an infinite series 
and includes correction terms for frequency (due to 
the finite time of propagation), as well as corrections 
for the finite radii of both the transmitting and receiv- 
ing loops. The derivation is based on a current of 
constant amplitude and [)hase around the transmitting 
loop. However, the effect of a nonuniform current dis- 
tribution at fre(piencies above 5 or 10 MHz is also 
analyzed. Formulas used by other workers in the field 
[15] to 1 17] often fail to include these refinements 
which can result in errors of up to 20 percent and more 
in computing standard-field values. 

The accuracy of the new fornuila is basically limited 
only by the number of terms used in the series. At the 
highest fre([uency normally emph^yed, 30 MHz, tlu^ 
error in one of the approximate forms given does not 
exceed 0.2 percent for the conditions stated and de- 
creases rapidly with decreasing fre(piency. \\w series 
expression is rapidly converging and simple to use by 
workers in the field without recourse to a table of 
functions or a computer. 
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Bessel functions of odd half order have particularly 
simple forms in terms of trignometric functions [22]. 
A few of these iclalionships are as follows: 



J\l2[x) = \j — sin X 



,2 / sin X 

J3/2(X) = \ COS .X- 

7TX \ X 



(A2) 



(A3) 



J^l'zix)-- 



2_ 

7TX 



sm X cos X I (A4) 

X 



Nij'zix) = \/— - (— cos x) 

TTX 



M ( \ / 2 / . COS 

I^'^l2{x) = \ — I — sm a: 

V TTX 



/ 2 / 3 
N:,j2{x)=\ — sin X- 

y TTX \ X 



)S x\ 
X / 

3__ 

x'- 



(A5) 



(A6) 



cosjf . (A7) 



The recurrence relations for calculating these func- 
tions for higher ord(Ms are [23 1 



7/1 + 1 {x)= — J„{x) -Jn-i{x) 

X 



2n 



(A8) 



Nn^Ax)=—Nn{x)-iyn-l{x). (A9) 

If these two relations are substituted into the definition 
of" the s[)herical Hankel functions, given l)y ecj (Al) of 
this appendix, the recurrence I'clation foi- spherical 
Hankel functions is obtained as [24] 

/^/ + i(-^) =— hn{x)-hn-i{x). (AlO) 

Using the above relations several spherical Hankel 
functions of higher order can be determined. 



8. Appendix 

For the convenience of the reader, explicit expres- 
sions for the various Bessel functions used herein, but 
not commonly found in compilations, are listed below. 

The spherical Hankel functions of the second kind, 
used in this paper, are defined in terms of Bessel func- 
tions of the first and second kinds, J„{x) and A^, (^) 
respectively, by the relation [21] 



(Al) 



h^~\x) - V ^ ^'1'^ + ^/^ W -y/V. + i/2(x)), 



h^^Kx)=^^(^j) 

X 



/7<2'(:r) = ^-^(-l+- 



X \ XX 



X \ XX- : 



(All) 
(A12) 
(A13) 

(A14) 
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I>?{x) 



hfHx) 



X 


\ . , 10 .45 105 , .105 

^ X X' X'^ x^ 


X 


' , , .15 , 105 .420 94 

^ X X~ X'^ x^ 



(A15) 



420_945^ . 945 \ 



(A16) 

It should be noted that the expression for each of 
the above functions is exact, each being a polynomial 
of degree (ai+1). These Hankel functions can also 
be derived somewhat more directly from the following 
relation based on Erdelyi [25]: 



A<,?'W=;(-l)"^"(4- 
\xax. 
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